Pseudoscalar-meson-Octet-baryon Coupling Constants in Two-flavor Lattice QCD 
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We evaluate the nNN, 7rE£, 7rAE, KAN and KT.N coupling constants and the corresponding 
monopole masses in lattice QCD with two flavors of dynamical quarks. The parameters representing 
the SU(3)-flavor symmetry are computed at the point where the three quark flavors are degenerate 
at the physical s-quark mass. In particular, we obtain a = F/(F + D) — 0.395(6). The quark-mass 
dependences of the coupling constants are obtained by changing the u- and the d-quark masses. We 
find that the SU(3)-flavor parameters have weak quark-mass dependence and thus the SU(3)-flavor 
symmetry is broken by only a few percent at each quark-mass point we consider. 
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Meson-baryon coupling constants are important ingre- 
dients for hadron physics as they provide a measure of 
baryon-baryon interactions in terms of One Boson Ex- 
change (OBE) models, and production of mesons off 
the baryons. In phenomenological potential models, the 
meson-baryon coupling constants are determined so as 
to reproduce the nucleon-nucleon, hyperon-nucleon and 
the hyperon-hyperon interactions in terms of, e.g., OBE 
models. On the other hand, it is an important issue to 
determine the coupling constants at the hadronic vertices 
directly from QCD, the underlying theory of the strong 
interactions. The only method we know that provides a 
first-principles calculation of hadronic phenomena is lat- 
tice QCD, which serves as a valuable tool to determine 
the hadron couplings in a model-independent way. 

Among other meson-baryon coupling constants, the 
nNN coupling constant, g^NNi which enters as a fun- 
damental quantity in low-energy dynamics of nucleon- 
nucleon and pion-nuclcon, has been a subject of intense 
investigation. It is defined as the nNN form factor, 
9-kNn{<I 2 ), at zero momentum transfer, q 2 = 0. The value 
of the coupling constant at the pion pole is relatively 
well-known from experiment: g\ NN (rri^) j Ai: ~ 13.6 (see, 
e.g., Rcf. pj, for a review). The value at zero mo- 
mentum transfer can be extracted from the Goldberger- 
Trcimann relation (GTR), g^NN = Qa^nI ' f-n ~ 12.8, 
where / ff is the pion decay constant and and qa are 
the mass and the axial-vector coupling constant of the 
nucleon, respectively. An earlier determination of the 
7riVJV coupling constant from a quenched-lattice QCD 
calculation, which reports g^NN = 12.7 ± 2.4 3, is 
in agreement with the phenomenological value. In the 
SU(3)-fiavor [SU(3)f] symmetric limit, one can classify 
the pseudoscalar-meson-octet-baryon coupling constants 
in terms of two parameters: the nNN coupling constant 
and the a = F/ (F+D) ratio of the pseudoscalar octet 
This systematic classification is expected to govern all the 
meson-baryon couplings, however as we move from the 
symmetric case to the realistic one, the SU(3) f breaking 
occurs as a result of the s-quark mass and the physical 
masses of the baryons and mesons. The broken sym- 



metry no longer provides a pattern for the meson-baryon 
coupling constants, and therefore they should be individ- 
ually calculated based on the underlying theory, QCD. 

In this work we extract the coupling constants g^NN, 
5irss, ftrAE, 9k an and g K SN (denoted by g M BB> here- 
after) by employing lattice QCD with two flavors of dy- 
namical quarks. The evaluation of the coupling constants 
allows us to check the following SU(3) f relations: 



9vNN — 9, 
9k an = 



1 



25" i 5?rA£ 



.9(1 - a), 



(1) 
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g(l + 2a), gxEjv = .9(1 - 2a), 



which phcnomenologically work rather well but are not 
known a priori to hold. 

The pseudoscalar current matrix element is written as 

(S(p)|P(0)|B'(p')) - gp(q 2 )u(p)i l5 u(p'), (2) 

where gp(q 2 ) is the pseudoscalar form factor, q^ = 
p' — p^ is the transferred four-momentum and P(x) = 
i/5(x)i75 : y'0(x) is the pseudoscalar current. We compute 
this matrix element using the ratio @] 



(GB'(t 2 ;p' ; r 4 )) 



(G B (t 2 -t 1 -p;T 4 )} 
(GB'(t 2 -t i; p';r 4 )) (3 ) 



(G g 'ft i; p^r 4 ))(G g '(t 2 ;p';r 4 )) 
(G^ i; p;r 4 ))<GS(t 2 ;p;r 4 )) 



1 1/2 



where the baryonic two- and three-point correlation func- 
tions are respectively defined as 



(G e (t;p;r 4 ))=^e-^rr' 

X 

x (vac|T[7 ? g( a; )77g;(0)]|vac), 
{G BVB \t 2 ,t 1 ;pi,p;T)) = - l £ 



3 -ipx 2 „iq-xi 



(4) 
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x T aa (vac|T[r;g(x 2 )P(x 1 )^,(0)]|vac) 
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FIG. 1: The q dependence of the form factors, Qmbb' for 
K v'ai = 0.1393. The diamonds show the lattice data, and the 
solid curves denote the fitted form factors. 



with r = 7375^ and r4 = (1 + 74)/2. The baryon 
interpolating fields are given as 



m (x) = e abc [u Ta (x)C 75 d b (x)]u c (x), 
Mx) = e abc [s Ta (x)C l5 u b (x)}u c (x), 

VA (x) = ^=e abc {[u Ta {x)C l5 s b {x)]d c {x) 



(G) 



[d Ta (x)C 



x l5S b (x)]u c (x) + 2[u Ta (x)C l5 d b (x)]s c (x)}, 
where C — 7472 and a, b, c are the color indices. t\ is 



FIG. 2: g-xNN and A^jvat as a function of m%. The empty 
circle denotes the SU(3)f limit and the diamond marks the 
experimental result. 
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FIG. 3: The ttSS, ttAS, KAN and KEN coupling constants 
normalized with g-^NN as a function of m%. The empty circle 
denotes the SU(3)f limit. 



the time when the meson interacts with a quark and ti is 
the time when the final baryon state is annihilated. The 
ratio in Eq. ([3]) reduces to the desired pseudoscalar form 
factor when ti — t\ and t\ 3> a, viz. 



i?(t 2 ,t i; 0,p;r;/i) 



ti»a 



g L P (q 2 ) 



■93, 



(7) 



t 2 -t 1 >a [2E(E + m)} 1 / 2 
where m and E are the mass and the energy of the ini 
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tial baryon, respectively, and gp(q 2 ) is the lattice pseu- 
doscalar form factor. Since the ratio in ([7]) is propor- 
tional to the transfered momentum q% . it cannot be used 
directly to obtain gji(q 2 ) at q 2 = 0. We apply a proce- 
dure (similarly to the one in Ref. [U) of seeking plateau 
regions as a function of t\ in the ratio ([7|) and calculating 
gp(q 2 ) at the momentum transfers q 2 a 2 = n(2ir/L) 2 (for 
the lowest nine n points), where L is the spatial extent 
of the lattice. We then obtain the meson-baryon form 
factor via the relation 



9 L P {q 2 ) = 



Gm Qmbb 1 (q 2 ) 



(8) 



M 



assuming that the pseudoscalar form factors are dom- 
inated by the pseudoscalar- meson poles. Here Gm = 
(vac|P(0)|M) is extracted from the two-point mesonic 
correlator (P(x)P(O)). Finally we extract the meson- 
baryon coupling constants guBB 1 = Qmbb 1 (0) by means 
of a monopolc form factor: 



gMBB'(q 2 ) = guBB 1 



A 2 

2 ^MBB> 



A 2 

MBB' 



(9) 



We employ a 16 3 x 32 lattice with two flavors of dynam- 
ical quarks and use the gauge configurations generated 
by the CP-PACS collaboration Q with the renormaliza- 
tion group improved gauge action and the mean-field im- 
proved clover quark action. We use the gauge configura- 
tions at (3 = 1.95 with the clover coefficient c$w = 1-530, 
which give a lattice spacing of a = 0.1555(17) fm (a -1 = 
1.267 GeV), which is determined from the p- meson mass. 
The simulations are carried out with four different hop- 
ping parameters for the sea and the u,d valence quarks, 
K sea , «w = 0.1375, 0.1390, 0.1400 and 0.1410, which cor- 
respond to quark masses of ~ 150, 100, 65, and 35 MeV, 
and we use 490, 680, 680 and 490 such gauge config- 
urations, respectively. The hopping parameter for the 
s valence quark is fixed to k. 



vol 



0.1393 so that the 
Kaon mass is reproduced Q, which corresponds to a 
quark mass of ~ 90 MeV. We employ smeared source 
and smeared sink, which are separated by 8 lattice units 
in the temporal direction. Source and sink operators are 
smeared in a gauge-invariant manner with the root mean 
square radius of 0.6 fm. All the statistical errors are es- 
timated via the jackknife analysis. 



TABLE I: The fitted values of m T , ma, tun, mh and mj in 
lattice units. 



niK 



TUN 



0.1375 0.899(1) 0.834(1) 1.707(06) 

0.1390 0.737(1) 0.725(1) 1.475(05) 

0.1393 0.713(1) 0.713(1) 1.455(06) 

0.1400 0.603(1) 0.635(1) 1.289(05) 

0.1410 0.440(1) 0.533(1) 1.051(08) 



1.658(06) 1.648(06) 

1.466(06) 1.464(06) 

1.455(06) 1.455(06) 

1.312(04) 1.318(05) 

1.114(06) 1.134(07) 



In Table U we give the fitted values of the meson and 
baryon masses as obtained from the two-point correlation 
function in Eq. (j4|). We extract the meson-baryon cou- 
pling constants, gMBB' , and the corresponding monopolc 
masses, Ambb', for each k"^. In Fig. [TJ the q 2 depen- 
dence of the form factors, guBB 1 ■ for K v' a i = 0.1393 
is given. Our complete results are presented in Ta- 
ble HU We give the fitted value of the irNN coupling 
constant and the corresponding monopole mass in lat- 
tice unit, as well as the fitted values of the 7rEE, 7rAS, 
KAN and K'SN coupling constants and the correspond- 
ing monopole masses normalized with g^NN and A^at, 
respectively. In Table [III 9aibb' an d A{f /BB , denote 
gMBB'/g^NN and Ambb'/A^nn, respectively. We ex- 
pect that the systematic errors cancel out to some degree 
in the ratios of the coupling constants and those of the 
monopolc masses. We give a graphical representation of 
our results in Figs. [21 El and [U In Fig.[5]we plot g^NN and 
A-nNN as a function of the pion-mass squared. The ratios 
of the 7rES, 7rA£, KAN and KUN coupling constants 
to the kNN coupling constant, and the corresponding 
monopolc masses normalized with A^nn are shown in 
Fig- El gnNN is consistent with the experimental value 
at k < 0.1393 and A^nn decreases towards the chiral 
limit. Note that, in addition to the monopolc form, we 
have tried fitting the form factors to dipole and exponen- 
tial forms, which have produced coupling-constant ratios 
consistent with those given in Table HH 
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FIG. 4: Same as Fig.[3]but for monopole masses A^es, A^as, 
Akan and Aksn normalized with A w nn- 
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TABLE II: The fitted value of the irNN coupling constant and the corresponding monopole mass (given in lattice units), 
together with the fitted values of the 7r££, 7rA£, KAN and KT,N coupling constants and the corresponding monopole masses 
normalized with g-^NN and A-n-Arjv, respectively. Here, we define g^ggi = Qmbb 1 / 'g^NN and A R fBB , — Ambb' / J^-ttNN ■ 



K val 


g-KNN 


A,rAfJV 


n R 
<?ttEE 


n R 


9kan 


9ksn 


\ R 


A ft 

a ttAE 


A JfAiV 


\ R 



0.1375 13.953(412) 1.053(123) 0.759(11) 0.698(11) -1.038(07) 0.231(14) 1.074(065) 0.908(039) 1.011(27) 0.714(118) 

0.1390 13.257(448) 1.228(189) 0.785(12) 0.697(07) -1.034(07) 0.209(12) 1.020(066) 0.988(042) 1.006(28) 0.978(223) 

0.1393 13.236(478) 1.223(202) 0.789(13) 0.699(08) -1.033(08) 0.209(13) 1.020(068) 0.989(044) 1.009(30) 0.970(236) 

0.1400 13.098(393) 1.013(111) 0.781(13) 0.723(08) -1.017(07) 0.242(15) 1.034(053) 0.970(033) 1.026(24) 0.802(124) 

0.1410 12.834(1.092) 0.719(123) 0.781(38) 0.756(28) -1.007(30) 0.260(30) 1.083(106) 0.985(074) 1.032(58) 0.958(191) 



Having discussed the results for g-^NNi we proceed 
with the octet-meson-baryon coupling constants. Wc 
first concentrate on the SU(3)-flavor symmetric case, 
where K™' al = K^ al = 0.1393 and the SU(3)^ relations 
in Eq.([T]) are exact. (Here we take k™^ =0.1390 and 
neglect the difference in the sea-quark effects.) As ex- 
pected, all the coupling ratios, e^ EE , ^ AE , g% AN , and 
9khn are wen reproduced with a = 0.395(6), which is 
obtained by a global fit. The ratios of the monopole 
masses, A^ ss , A^ AE , A^ AAr , and A KSN , are consistent 
with unity. The obtained value of a is consistent with 
that in the SU(6) spin-flavor symmetry (a = 2/5) [H, 
which is the symmetry based on the nonrelativistic quark 
model. We have also tried fixing all the quark masses 
at k"^' s = 0.1390. g-wNN and the ratios of the cou- 
pling constants obtained in this case are as follows: 
g nNN = 12.769(495), <?« ES = 0.785(10), ff « AE = 0.704(6), 
9k an = -1-003(6), and g§^ N = 0.211(10). We have 
found that the coupling constants again satisfy SU(3)i? 
and the resulting a = 0.387(5) is consistent with that 
obtained at k"^' s = 0.1393. In Fig. [5] we plot the ra- 
tio in Eq. §3§ for g^NN as a function of current-insertion 
point, in order to show the plateau regions. We present 
the data at «"^' s = 0.1393 for the first three momentum 
transfer values. 
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FIG. 5: The ratio in Eq. (|3]) for g W NN as a function of current- 
insertion point, ti, at k"^' s = 0.1393 for the first three 
momentum-transfer values. The bands represent the adopted 
plateau regions. 
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We next discuss the SU(3) p broken case. The quark- 
mass dependences we find for g^jBB' ano - ^mbb' are n °t 
large. The ratios of the coupling constants, g^BB 1 ' arc 
similar in value to those in the SU(3)f symmetric limit, 
and the monopole-mass ratios, A^ [B b' > are almost unity 
independently of the quark masses. This suggests that 
the SU(3)f breaking is small at the quark masses we 
consider. Our data do not allow a direct determination 
of SU(3)f breaking in the chiral limit, as wc have flavor- 
symmetric data only at K%' al = n s val = 0.1393. On the 
other hand, the value of g n NN at the chiral point is well- 
known, which may serve as a reference point for us to 
obtain a measure of SU(3)f breaking. For this purpose, 



FIG. 6: The value of 5su(3) as a function of m\. The curve 
and the shaded region denote linear chiral extrapolations with 
errors. 

we construct the following three sets of relations: 
Al =2 (^ 9 " AS + ' M = 9 ^ N + 

A 3=2 { 9 ^ N ~ ^9kAn) > A i = ~g%T, - V^fff AJV , 

A =-= (g* AS - g KKN ) , A 6 = Vdg^ - g% SN , (10) 
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FIG. 7: The value of a as obtained from a global fit of the 
SU(3)f relations at each quark mass we consider (in black 
filled circles). We also show our results as obtained with local 
source and local sink (in red triangles). The empty circle and 
the triangle denote the SU(3)f limit. The line at a = 0.4 is 
shown for reference only. 



and 



B 2 
B 3 

Bi 



AS + 3ff^ss 



^9kt,n 



~ ' + ^9ksn 



n R 



4.As 



^9kan 



^9kan^ 
.9^ss^ 
9ky,n — 9k AN I ! 



V3q R 



(11) 



V3r R 



5ttAS 



V3r R 



9k AN 



n R 

9tt ss 



9ksn 



(12) 



which can be readily obtained from those in Eq. JTJ). In 
the SU(3)f symmetric limit, the above equations sat- 
isfy A x = ... = M = B x = ... = B A = d = 1, 
which can be verified by inserting the coupling con- 
stants at K™' al = 0.1393 in Table [TTJ At other quark 
masses, the deviations from unity represent the amount 
of SU(3) f breaking. Inserting the values of the coupling 
constants corresponding to the lowest quark mass we con- 
sider in Table OH into (Hni)-(II2>, we find A 1 = 1.045(29), 
A 2 = 1.040(30), A 3 = 1.002(25), A A = 0.963(42), 
A 5 = 1.017(22), A 6 = 1.049(40), B x = 1.043(28), 
B 2 = 1.021(24), B 3 = 0.990(28), B 4 = 1.033(28) and 
C\ = 1.006(27), which indicate a breaking in SU(3)f by 
less than 10%. Moreover, we define the average SU(3)f 
breaking as follows: 



3 SU(3) 



1 

IT 



E i 

n,X=A.B,C 



i-x n \, 



(13) 



which amounts to <5 SU(3) =0.014(03), 0.003(02), 
0.012(02), and 0.028(17) for the quark masses at - 150, 
100, 65, and 35 MeV, respectively. This suggests for the 
pscudoscalar-meson couplings of the octet baryons that 
SU(3) f is a good symmetry in the quark-mass range we 
consider, which is broken by only a few percent. We 



have also tried a quadratic fit of ^su(3) and extracted 
<5su(3) = 0.072(16) in the chiral limit. Fig. [6] shows the 
value of <5gu(3) as a function of and the chiral extrap- 
olations with errors. In Fig. [71 we plot the value of a as 
obtained from a global fit of the SU(3)i? relations at each 
quark mass we consider, a slightly decreases toward the 
chiral limit . The deviation of a in the present quark-mass 
range is at most 10%, whereas that of c5gu(3) is less than 
5%. We infer from this that the deviation in a should 
be small in the chiral limit, as we find that the ratios 
of the coupling constants have weak quark-mass depen- 
dence. The SU(3) f breaking effect seems to appear in a 
rather than in SU(3)f relations (£su(3))- We have also 
repeated our analysis with local source and local sink for 
consistency check. In Fig. [JJ we show the values of a as 
obtained from such a setup as well, where both analysis 
lead to consistent results with each other. 

In summary, we have evaluated the pseudoscalar- 
meson-octet-baryon coupling constants, g^NN, ff-n-ss, 
.9ttAS, 9 kan and <?a'sw, in two-flavor lattice QCD with 
the hopping parameters n sea , n™' al = 0.1375, 0.1390, 
0.1400 and 0.1410, which correspond to quark masses of 
~ 150, 100, 65, and 35 MeV. The parameters represent- 
ing the SU(3)f symmetry have been computed at the 
point where the three flavors are degenerate at the phys- 
ical strange-quark mass. In particular, we have obtained 
a = F/(F + D) = 0.395(6), which is consistent with the 
prediction from SU(6) spin-flavor symmetry (a = 2/5). 
The monopole mass we find leads to a 7riV/V form fac- 
tor which is softer than those typically used in the phe- 
nomenological OBE potential models. The ratios of the 
coupling constants, which are supposed to be less prone 
to systematic errors, show very weak quark-mass depen- 
dence. We have discussed to what extent the SU(3)f 
symmetry is broken as we approach the physical masses 
of the u- and the (i-quarks. Our results indicate for the 
pseudoscalar-meson couplings of the octet baryons that 
SU(3) f is a good symmetry, which is broken by only a 
few percent (at least) in the 35 MeV to 150 MeV range 
of the light quark masses. 
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